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AN ANALYSIS OF THE RADIATION FIELD BENEATH A BANK 
OF TUBULAR QUARTZ LAMPS 
By Robert L. Ash^ 

SUMMARY 

Equations governing the incident heat flux distribution beneath a lamp- 
reflector system have been developed. Analysis of a particular radiant heating 
facility showed good agreement between theory and experiment when a lamp power 
loss correction was used. In addition, the theory was employed to estimate 
thermal disruption in the radiation field caused by a protruding probe. 

INTRODUCTION 

One method for producing relatively high surface heating rates has been 
to employ banks of tubular, tungsten filament quartz lamps. These lamps are 
commercially available in lengths^ of between 12.7 cm and 96.5 cm and are capa- 
ble of dissipating up to 35 w/cm (ref. 1). By employing a group of these lamps 
beneath a reflecting panel, surface heating rates in excess of 50 w/cm^ can be 
developed. 

The purpose of this report is to present a theory which can be used to 
estimate the incident energy distribution produced by lamp-reflector systems. 

The configuration under study is shown schematically in figure 1 (p. 2). 

For simplicity, it was assumed that the test panel is perfectly absorbing 
(black) and is at a sufficiently low absolute temperature to neglect radiation 
back to the reflector. Furthermore, it was assumed that the lamp diameters are 
small relative to other dimensions of the system and that they do not interfere 
with the radiant energy passing between reflector and test panel. 


1 Associate Professor of Engineering. School of Engineering, Old Dominion 

University, Norfolk, Virginia 23508. ^ , a.- * n -jo's? Anrh 

2 Length measurements appear in the metric system, 1 centimeter - 0.3937 men 



Figure 1.- Sahematic view of radiant heating system. 


ANALYSIS 

It is convenient to focus attention on a single lamp having an arbitrary 
location beneath a flat reflector panel. Analysis of the radiant energy leaving 
the lamp can be divided into two parts, direct radiation of the test panel and 
reflected radiation to the test panel. Direct radiation calculations will be 
developed first since they can be used to simplify the reflector analysis. 

Because the lamp radiates uniformly about its axis, it should be realized 
that the unit normal on a differential cylindrical surface element is free to 
rotate in the plane shown in figure 2 (p. 3). A complicated integration can be 
avoided by using the indirect approach of assuming that the intensity leaving 
the cylindrical surface element of figure 2 is some known function Ii(y). 

That is, Ii(y) given by 

ii(y) = 

0 

is treated as a known quantity, where Il(T) is the Intensity leaving a dif- 
ferential surface element (rectangular) of the lamp at temperature T and 
f( 0 ) is the function which relates the intensity leaving the differential 
element to the particular direction under consideration (f(0) will be zero for 
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Figure 2,- Plane of radiation eynvnetry. 


some angles since part of the surface is hidden by the lamp filament). It will 
be convenient to evaluate IjCy) later; but it should be noted that, if the 
lamp is at a uniform temperature, Ii(y) would be a constant— say 1^. 

There are several advantages to this indirect approach. First of all, 
as has already been mentioned, a difficult integration will at least be post- 
poned. More importantly, by taking advantage of the symmetry of the radiant 
energy leaving the lamp, it will be possible to Ignore the radiation field 
within the quartz cylinder generated by the tungsten filament. No transmittance 
calculations are needed. In fact, when the uniform filament temperature assump- 
tion is made, it will be found that the filament temperature itself is not 
needed in the calculations. More precisely, it will be possible to relate 
directly I^ to the total power dissipated by the lamp. 




where D is the lamp diameter. This equation can be rewritten in terms of the 
dimensionless variables f 2 , I, and y» defined by 
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The virtual intensity can be related to the lamp power by realizing that 
one-eighth of the total energy radiated by the lamp (Pq/S) is incident on the 
test surface region given by 0 ^ y < «>, 0 £ x < ®. 

Hence, 
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Equation (5) can be integrated directly to yield 
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Now. when the lamp is at a uniform temperature and there is no radiation from 
the surroundings. 
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Furthermore, 


frl = eoT'^ 
0 0 


( 10 ) 


Though these last two relations are interesting, their usefulness is 
limited since they relate temperature to heating rate and yield an apparent 
emissivity for the cylindrical element when the temperature is known. There 
is no well defined temperature or emissivity because both quartz cylinder and 
tungsten filament emit radiant energy at different temperatures with different 
emissivities, to say nothing of the transmittance through quartz. 

If the nth lamp is located at a distance from the reference edge 


of the reflector, x is related to the reflector coordinate x by 
X = X - x„ 


n 


Similarly, 


where 


and 


C = C - 5. 
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( 11 ) 
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the incident flux produced by the nth lamp is given by 


fn(^*y) = ^ F(x,y;x^,H) 


where 
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, and ? are defined by 
iin = [(x - x^)2 + H^JVh 


Y = L/H 


and 


c = y/L 

F(x,y;x^,H) shall be called the radiation function. 


(13) 
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Reflected Radiation Contribution 

Reflected contribution to the energy distribution on the test panel can 
be developed conveniently from the previous calculations. If the reflector is 
assumed to act as a combination specular and diffuse surface, the reflec- 
tivity p can be written 

(17) 


p = 


d . s 

p + p 


where p^ is the specular reflection component and p'^ is the diffuse compo- 
nent. Separate analyses are required for the specular and diffuse 
contributions. 

The specular contribution can be developed directly from the previous 
direct radiation analysis. That is, if the distance from the lamps to the re- 
flector panel is R, a set of imaginary lamps located at a distance H + 2R 
above the test panel and having strength given by 
s- 


's = 


pi 


(18) 


represent the radiant energy caused by specular reflection from the test panel. 

If the specular contribution of the nth lamp to the incident flux 
at the test panel is designated f^(x,y), it should be obvious that 


f"(x,y) = 




, -F(x,y;x ,H + 2R) 

(H + 2R)2 ^ 


(19) 


In order to calculate the diffusely reflected energy contribution, it is 
first necessary to calculate the incident energy distribution on the reflector 
panel. Again, the incident energy function can be developed directly from the 
radiation function. In this case, the teat surface is located at a distance R 
from the lamps and the incident flux on the reflector panel is given by 


Bn(x»y) = pF(x,y;x^,R) 


( 20 ) 


Consequently, the intensity leaving a point on the reflector panel would be 
p‘^Bj^(x,y)/n and the incident flux on the test panel due to diffuse reflection 
from the nth lamp would be 
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where 0|^j and are the angles between normals at the reflector and test 

panel and the radius r^j^. From this particular geometry, 


cos = cos 0-pp = (R + H)/r 


where 


TR 


r|p = (R + H)2 + (x - u)2 + (y - v)^ 


( 22 ) 


(23) 


From equations (22) and (23), the incident energy on the test surface can be 
wri tten 

d 
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Because of the form of F(u,v;x^,R) in equation (15), equation (24) is 
a hyperel 1 iptic integral and must be integrated numerically. A computer pro- 
gram has been written to evaluate equation (24), using an equally spaced grid 
network. The integral has been approximated by 
(R . h) 4 f F(u.v;x„,R)dud v 

^ JJ[(R + H) + (x - u)2 + 


F^(x,y;x^,R,H) 
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(25) 


i=t J=l[(R + H)2 + (x - U.)2 + (y - Vj)2]2 

where F^ has been defined for convenience and Nj and Nj represent the 
number of reflector surface elements in the x and y directions, 
respectively. 

At this point, the incident heat flux due to the nth lamp located at 

some distance x from the reflector reference line can be written 
n 
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Finally, the incident heat flux caused by combined radiation from an array of 
N lamps is given by 

q(x»y) = I fl(x.y) 

n=l 


N 

— y 

^ n=l 


+ YjP 


|^F(x,y;x^,H) + Y5p*F(x,y;x^,H + 2R) 
"^Fj,(x,y;x^.R,H)1 


(27) 


An experimental investigation has been carried out to determine the ac- 
curacy of this model and will be discussed further in the experimental section. 
A computer program for evaluating equation (27) when p = 0 is given in 
appendix A. 


EXPERIMENT 

The lamp-reflector system studied was located at the National Aeronau- 
tics and Space Administration, Langley Research Center and was used for the 
study of space shuttle material candidates. The system was designed to produce 
a nearly uniform radiant heat flux over a 5-cm by 5-cm area in an effort to 
simulate multiple reentry conditions. 

Essentially, the system consisted of 45 cylindrical, tungsten filament 
lamps mounted beneath a water-cooled aluminum reflector similar to the sche- 
matic shown in figure 1 (p. 2). The lamp-reflector system was mounted on a 
water-cooled aluminum table and could be adjusted vertically on the table. The 
lamps were mounted in water-cooled copper brackets and equally spaced 1.27 cm 
apart. The distance between the plane of the lamp centerlines and the reflector 
was 3.3 cm. 

Power was supplied to the lamps by a balanced three-phase system which 
supplied each group of 15 lamps (1-15, 16-30, and 31-45) with nominally the 
same power. The power consumption was measured by monitoring the phase voltage 
and current. 

The incident heat flux was measured with calibrated thin-foil, Gardon- 
type, heat flux sensors with foil diameters of 0.635 cm and surface emissivities 
of 0.89. The sensors were mounted in a low thermal -conductivity block, which 
positioned the sensing foil 5.08 cm above the water-cooled table. The heat flux 
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sensors were positioned on the horizontal plane by means of a grid system which 
had been inscribed on the surface of the table. 


The major difficulty in comparing theory with experimental data rested in 
separating the energy which was radiated by the lamps from the energy con- 
ducted through the lamp mounts and from the energy carried away convectively 
by the air. Only the radiant energy leaving the lamps was included in the 
theory. 

Initially, energy balances were made on the cooling water passing through 
the copper lamp mounts to determine conduction loss. However, cooling water 
data were not accurate because the brackets were heated by radiation and were 
receiving energy by conduction from the lamps. The final procedure consisted 
of measuring the conduction-convection loss experimentally. 


The experimental procedure employed assumed that theory and data agreed 

at a selected data point. By forcing agreement at the particular data point, 

a (power radiated by each lamp) could be calculated and the difference 

between the calculated and the power supplied to each lamp 

/total power supplied defined as the unradiated energy or power loss. 

\ total number of lamps/ 

Unfortunately, power loss was affected by the spacing between the lamp- 
reflector system and the water-cooled table. As the distance between the 
table and lamps was decreased, the lamp temperatures increased and the conduc- 
tion and convection losses increased. In addition, radiation loss resulting 
from multiple reflections between table, which was not perfectly absorbing, and 
reflector changed with spacing. All of these factors affected the calibration. 


Since most of the tests conducted at this facility positioned the speci- 
men surface either 8.9 cm or 5.08 cm from the test surface, conduction cor- 
rections were developed for these distances. The data were gathered when the 
heat flux sensor was located at 8.9 cm from the lamp plane and a detailed 
check was made at a single power setting at 5.08 cm. The conduction correction 
data are shown in figure 4 (p. 11). 
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Figiive 4.- Pouev loss oorreation. 


In order to verify the validity of the conduction correction, additional 
measurements were taken with the heat flux sensor in the reference position, 
but varying the lamp configuration. Using a power setting which nominally 
supplied 495 total watts to each lamp when the lamps were 5.08 cm above the 
sensor surface, two tests were run with 30- and 33-lamp configurations in place 
of the original 45. The 30- lamp configuration was produced by removing 15 lamps 
from their brackets so that there were 10 lamps equally spaced 1.27 cm apart, 
a gap of 7.62 cm, 5 lamps 1.27 cm apart, a gap of 7.62 cm, 5 lamps 1.27 cm 
apart, a gap of 7.62 cm, and 10 lamps 1.27 cm apart. The 33-lamp configuration 
was produced by replacing one lamp at the center of each of the 7.62-cm gap. 

The three test configurations are shown schematically in figure 5 (p. 12) 
along with the theoretically predicted heat flux distributions. 


11 




HORIZONTAL DISTANCE, cm 

Figure 5.- Beat flux measurements and distributions for calibration experiments. 

The data point marked reference in figure 5 was used to develop the power 
loss correction shown in figure 4 (p. 11). The calculated 45-lamp power loss 
was 218 watts per lamp. Hence, the value of for each of the lamps was 
277 watts. Using that value in the theoretical calculations, distributions 
shown in figure 5 were generated. The differences between predicted and meas- 
ured heat fluxes in the 30- and 33-lamp configurations were -1 percent and 
+2.5 percent, respectively. 

Obviously, the above correction procedure makes the theory appear more 
accurate than it is in reality. However, the data gathered in a previous test, 
which employed the 45-lamp configuration 5.08 cm above the test surface, verify 
the fact that measured and predicted incident flux distributions are in close 
agreement when the power loss correction is used. The flux was measured in the 
X and y directions, 2.54 cm off the centerline in both cases. The power loss 
was taken directly from figure 4 and no attempt was made to match data. The 
measured and theoretical distributions are shown in figures 6 and 7 (p. 13). 

It can be seen that theory predicts both heat flux magnitude and distribution. 
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Figure 6.- Measured and theoretical heat flux distributions for a 45-lcanp test 
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Figure 7.- Measured and theoretical heat flux 
distributions for a 45-lamp test. 

(H =. 5.08 cm. X to X „ - 2.54 cm) 



From the experimental results discussed in this section, it can be con- 
cluded that theory and experiment show good agreement. The correction factor 
employed in the analysis probably overcorrects the theory and, therefore, re- 
stricts quantitative accuracy statements. However, data gathered for the 
three configurations indicated that the correction is valid, at least for 
these particular lamp configurations. 

THERMAL DISRUPTION DUE TO A PROTRUDING PROBE 

One proposed method for measuring surface temperatures in the previously 
described radiation facility utilizes a protruding cylindrical probe. The 
probe configuration is shown schematically in figure 8. 

PROSE 
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Figure 8.- Full-scale view of probe protruding through center of reflector. 

Essentially, a hole will be cut in the center of the reflector to allow 
the probe to extend through the radiation field. The probe will be lowered 
through the hole until it contacts the test surface and remain there long 



enough to record the surface temperature. Then, the probe will be removed to 
permit uniform radiant heating. The question is, how will the radiation field 
be affected in the vicinity of the probe during the time the probe is lowered 

in the field? 

When the probe is in contact with the test surface, no radiation is re- 
ceived directly beneath the probe. In addition, the surface in the vicinity of 
the probe is obscured from some parts of the lamp-reflector system. If the 
probe is assumed black and nonradiating, determination of thermal disruption is 
straightforward, though tedius. 

The previous program for calculating incident heat flux distribution has 
been modified to determine disruption caused by the probe. Essentially, the 
image of the probe, as seen from the center of a particular test surface ele- 
ment (an arbitrary surface element from the original computations), was projected 
back into the lamp and reflector planes. Location and dimensions of the ob- 
scured lamp segments as well as locations of the obscured reflector surface ele- 
ments were then computed.^ Then, the radiant energy arriving at the test 
surface element from the obscured lamp and reflector elements was computed and 
subtracted from the original undisturbed heat flux value. The complete program 
for making these calculations is included in appendix A. 

The central region of the particular configuration studied is shown ap- 
proximately to scale in figure 8 (p. 14). The lamps were 5.080 cm from the test 
surface. The probe had a major diameter of 0.660 cm, a barv'el length of 

0.508 cm, and a barrel diameter of 0,254 cm. The length of the taper from the 

major diameter to the barrel was 1.016 cm. In the computer calculations, the 
center lamp was removed from each of the three phases, leaving 42 lamps. The 
center lamp was removed from the center phase to permit the probe to be located 
in the center of the radiation field, while the center lamps were removed from 

the other two phases to balance the power consumption. The ratio of the 

disrupted heat flux to the undisturbed heat flux in the vicinity of the probe 
is shown in figures 9 and 10 (p. 16). 

Figures 9 and 10 indicate that significant variations in tlie incident 
heat flux are confined to a very small region in the vicinity of the probe. 

The radiation field is depressed only slightly at distances greater than about 
1 cm from the probe. Furthermore, the actual heat flux depression should be 

^ A reflector element was considered obscured if its center were obscured. 
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Figure 9,- Ratio of dis- 
turbed incident 
heat flu:;: to 
undisturbed, 
heat flux 
along the 
X axis in the 
vicinity of a 
probe whose 
major diameter 
is 0.66 cm. 

The probe is 
located at the 
center of the 
lamp-reflector 
system. 



Figure 10.- Ratio* of dis- 
turbed incident 
heat flux to 
undisturbed 
heat flux 
along the 
y axis in the 
vicinity of a 
probe whose 
major diameter 
is 0.66 am. 

The probe is 
located, at the 
center of the 
lamp-reflector 
system. 
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less than that shown because some radiant energy will be reflected from the 
probe to the test surface and that reflected radiant energy has not been 
included in the present analysis. 


CONCLUSIONS 

From comparison of experimental and theoretical results, it can be con- 
cluded that the present theory is a reasonable approximation of the actual ra- 
diation field produced by a bank of tubular lamps. The theory accurately 
predicted the dimensionless variation in heat flux with location for the lamp- 
reflector systems studied. However, a lamp power correction was required to 
yield close numerical agreement. It has been shown that this procedure may be 
useful in designing specified radiation heat flux distributions (appendix B), 
although additional work is required to develop two-dimensional design 
capabil ities. 

Finally, the theory has been used to show that heat flux distribution is 
only slightly altered in the vicinity of a protruding cylindrical probe. 
Accuracy of these calculations must be validated by experiment. 
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APPENDIX A 


COMPUTER PROGRAM 


The FORTRAN IV computer program for calculating incident heat flux from a 
lamp-reflector system is incorporated in this appendix. Included is the pro- 
gram for calculating reduced incident heat flux caused by a probe protruding 
through the center of the lamp-reflector system. 


//ODASH JOB (0977»0404*65»3)*'ROBERT'tMSGLEVEL-(2fO»»CLASS-T 
// EXEC FORTGCLG 
//FORTtSYSIN DD * 

REAL L , ^ 

DIMENSION DIR(92»44) ♦DDIF(92»AA) tXI (45) tHnA) 

C N»NUMBER OF LAMPS 

N»42 

C R. height OF REFLECTOR ABOVE LAMPS 


c 

c 

c 


c 

c 


c 


R»1.3 

L-LENGTH OF LAMP HEATING ELEMENT 


L“ 10» 

RHO=REFLECTOR REFLECTIVITY 
RHO=0 #95 

READ IN SPATIAL DISTRIBUTION OF LAMPS 
DO 2 I»liN 
READ<5*1) XI(I) 

1 FORMAT(F15*8) 

2 CONTINUE 

H»HEIGHT OF LAMPS ABOVE TEST SURFACE 


H * 2 • 

INITIALIZATION OF DIR AND DDIF 


DO 4 I»l»92 
DO 3 J»l*44 
D1R( I • J) 

DDIF< I # J) =0* 

3 CONTINUE 

4 CONTINUE 
WRITE(6»500)H»N 

500 F0RMAT(5X» 'H= ' »F6« 2 ♦ lOX » • N- SIS#//) 

CALCULATION OF THE DIRECT INCIDENT INTENSITIES ON THE TEST SURFACE 


AND THE REFLECTOR 

DX=GRID SPACING IN THE X-DIRECTION 
DX=.25 

DY»GRIO SPACING IN THE Y-DIRECTION 
DY=.25 


0X2*DX/H 


DXR=DX/R 

DY2*DY/L 

GAM-L/H 

GAR*L/R 

GAD=(H/R)**2 

DO 7 K«1»N 
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APPENDIX A - Continued 


C 

C 


c 


c 

c 


c 

c 


00 6 I»l*92 
DO 5 


CI = I 
CJ»J 

CI*Cl-.5 


CJ»CJ-.5 

X«CI*0X2 

Y»CJ*DY2-*55 

GA-GAM 

XL»XHK)/H 

call FICT(X*Y*XL*GA*On 

DIR( I »J) »DIR( I *J)+DI tT AC 

INCIDENT reflector INTENSITY IS CALCULATED BY TREATING IT AS 

TEST SURFACE 


X«CI*DXR 


GA»GAR 


XL=»XI <X) /R 

CALL FICT(X»Y»XL*GA»DI ) 

DDIF ( I I J ) =DDI F( I ♦ J ) ■♦'DI 

5 CONTINUE 

6 CONTINUE 

7 CONTINUE 

CALCULATION OF IRRADIATION AT SELECTED POINTS 
DO 11 I«ltA6 
DO 10 J*2*22*4 


A 


CI*I 

CJ«J 

CI*CI-t5 

CJ-CJ-.5 

X=CI*DX2 

X*AND^Y^ARE^THE DIMENSIONLESS COORDINATES OF THE POINT ON THE TEST 
PANEL 


DO 9 II»1»92 

DO 8 

CII=II 

CJJ«JJ 

CII=CII-.5 

CJJ=CJJ-*5 

XX=CI I*DX? 

XX AND YY ARE THE DUMMY VARIABLES USED IN THE REFLECTOR 
INTEGRATION ROUTINE 


PI =3.1A15925 
RHM=(1.+R/H)**2 
XDS=IX-XX)**2 
YDS=(GAM*(Y-YY) )**2 
PREF3GAD*GAM*RHM#RH0»DX2*DY2/PI 
DEN=(RHM + XDS-^YDS>**2 
REFL = PREF*D0IF( I I »JJ) /DEN 
DIR( I tJ)»DIR< I *J)+REFL 
8 CONTINUE 
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APPENDIX A - Continued 


9 CONTINUE 

10 CONTINUE 

WRITE(6tlOO)DlR( 1*2) «DIR(I *6) *01R(1 *10) tOIRil »14) fDIR( 1 *18) »DIR( I f 
122 ) 

300 FORMAT <E16*8) 

100 FORMAT <6E16*8) 

200 CONTINUE 

1 1 CONTINUE 
110 CONTINUE 

DO 1100 1»A6*53 
DO 1000 J*16*22 
IF( J-18)949»948»9A7 

947 IF( J-22)949*948»949 

948 IF( I-46)950»950»949 

949 CONTINUE 
CI-I 
CJ-J 

CI-CI-.5 
CJ-CJ-.5 
X»CI»DX2 
Y»CJ*DY2-i55 
DO 900 II»1*92 
DO 800 JJ=1*44 
CII=II 
CJJ=JJ 
Cl I*CI I-.5 
CJJ=CJJ-.5 
XX-CII*DX2 
YY»CJJ*DY2-i55 
XDS=(X-XX )**2 
YDS=<GAM*< Y-YY) )**2 
PREF»GAD*QAM*RHM#RHO»DX2*DY2/PI 
DEN= ( RHM+XDS+YDS )**2 
REFL«PREF*DDIF( 1 1 *JJ)/OEN 
, DIRtl »J)*DIR(1»J)-^REFL 

800 CONTINUE 
900 CONTINUE 

950 CONTINUE 
1000 CONTINUE 
1100 CONTINUE 

PROBE SHADOW PROBLEM 


PROBE DIMENSIONS 

A>^LcNGTH OF CONICAL SECTION BETWEEN MAJOR DIAMETER AND BARREL 
A ® 0 • 4 

B=BARREL LENGTH LENGTH (IN.) 

B»0.2 

C«BARREL DIAMETER (IN.) 

C*0.1 

D«MAJOR DIAMETER (IN.) 

PROGRAM IS OPERABLE ONLY IF D IS LESS THAN 0.50 INCH. 

0*0.26 


( IN. 
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C REDUCED INCIDENT HEAT FLUX WILL BE CALCULATED FOR I-A6-52* J-16-2 

RC»C/2# 

RO*D/2* 

HIN = H 

XCL»11*375 

DO 103 I«46»52 

DO 104 J«16t22 

IC0NT»0 

IWARN-0 

H»HIN 

CI»I 

CU»J 

X»(CI-0*5 »*0.25 
Y-(CJ-0.5 )*0.25-5.5 
WRITE(6.600)X.Y»DIR I I *J» 

600 F0RMAT(/#5X**X» ‘ . F 10 #4 * 5X * • Y- • .F10.4 . 5X ♦ *0R IG I NAL Q- * *E16»6) 
IF( 1-46) 105»106*105 
106 IWARN-1 

C IF PROBE DIAMETER EXCEEDS 0.50 INCH. LOGIC MUST BE CHANGED HERE. 

lF(J-22) 105.368.105 

368 DIR(I.J)*>0. 

GO TO 36 
105 CONTINUE 

IF(IWARN-l) 50. 51. 51 

51 DB»(Y«»2-RC*»2)**.5 
DA-(Y»#2-RD*»2 )*#.5 
THXY»0. 

GO TO 52 

50 DXYS= (X-XCL)**2-*-Y**2 
DB=(DXYS-RC**2)**.5 
DA=(DXYS-RD*»2>*#.5 
XMOYM=(XCL-X)/Y 
ThXY=ATAN(XMOYM) 

52 CONTINUE 
APB=(A+B)/DA 
THA=ATAN( APB) 

BP=B/DB 

THB=ATAN(BP) 

D0TA»D/l 2.*DA) 

C0TA=C/(2.*DB) 

THPA*ATAN(DOTA) 

THPB=ATAN(C0TA) 

THXAP=THXY+THPA 

THXBP=THXY+THPB 

24 RA=H/SIN ( THA) 

RB=H/SIN( THB) 

XAP=RA*C0S( THA) / ( ( 1 .+( TAN ( THXAP ) )**2 ) **.5 ) *TAN( THXAP ) 

XBP = RB*COS< THB) / ( ( 1 .-f I TAN ( THXBP ) )**2 )**.5 > *TAN( THXBP ) 

YAP=XAP/ ( TANI THXAP ) ) 

YBP=XBP/«TAN(THXBP) ) 

XAL=X-XAP 

XBL=X-XBP 

YAL»Y+YAP 
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250 


260 

25 


X09 

107 

111 

112 

13 

l^ 

160 

16 

17 

18 

19 

20 

22 

21 


162 

161 

15 

108 


122 

23 


YBL=Y+Y0P 

IFUWARN-l) 250*260»260 

THXAM=THXY-THPA 

THXBM=THXY-THPB 

XAM-RA*COS(THA)/( (l.-MTAN(THXAMJ 
XBM»R8»COStTHB)/M IfHTAN(THXBM) 

YAM»XAM/(TANtTHXAM) ) 
YBM=XBM/(TAN(THX0M) ) 

XAR-X-XAM 


)»#2)**.5)#TAN(THXAM) 
)##2 )*»*5)*TAN(THXBM) 


XBR>X-XBM 


YAR=Y+YAM 

YBR=Y+YBM 

IF( ICONT-l)25»26*26 
CONTINUE 

CALCULATION OF OBSTRUCTED LAMP COORDINATES 


NP*N/2 

DO 108 ILA»1#NP 
XLA=XI ( ILA) 

IF(XLA-XBL) 109 1 109 » 107 
Y1»YBL-MXBL-XLA)/ITAN(THXBP) ) 

GO TO 13 

IFIXLA-XAL) 111»111*112 

Yl=YAL+( XAL“XLA)#( YBL-YAL> /(XAL-XBL) 

GO TO 13 

Y1«Y+ { X-XLA) / ( TAN< THXAP) ) 

CONTINUE 

IF ( Yl-5* > lAtlA* 15 
CONTINUE 

IF(IWARN-l) 160i22»22 
IF(XLA-XBR) 16»16»17 
Y2=YBR+(XBR-XLA) /{TAN(THXBM) ) 

GO TO 20 

IF(XLA-XAR)18»18»19 

Y2 »YAR+ ( XAR-XLA ) # I YBR-YAR ) / < XAR-XBR ) 

GO TO 20 

Y2»Y+ I X-XLA) / ( TAN( THXAMIj N 

CONTINUE 

1F( Y2-5i ) 21i21*22 
Y2»5. 

CONTINUE 

CALL FAST ( X*Y.XLA»Y1*Y2*H#DLP) 

FUN*DLP 

IF(IWARN-l) 161»162»162 
DLP=2#*FUN 

DIRII*J)»DIR(I»J)-DLP 

CONTINUE 

CONTINUE 

ICONT=ICONT+l 

IFMCONT-l )122*122*23 

CALCULATION OF OBSTRUCTED REFLECTOR COORDINATES 
H»H+R 


GO TO 24 
CONTINUE 



APPENDIX A - Continued 


26 DO 26 II«1»46 

ci»n 

XRI«(CI-.5)*»25 

IF(XRI-XBU30«30*31 

30 YLl«YBL+(XBL-XRn/(TAN(THX0P) ) 

GO TO 34 

31 IF(XRI-XAU32*32.33 

32 YLl-YALMXAL-XRI ) # ( YBL-YAU / ( XAL-XBL) 

GO TO 34 

33 YLl-Y+IX-XRn/dANiTHXAPn 

34 CONTINUE 
IF(YLl-5.5>35. 35*36 

35 CONTINUE 
IF(IWARN-l) 370*43*43 

370 IF(XRI-XBR) 37*37*38 

37 YR2»YBR+( XBR-XRI)/ITAN(THXBM ) ) 

GO TO 41 

38 IF(XRI-XAR)39*39*40 

39 yR 2 = YAR-MXAR-XRI ) * t YBR-YAR ) / ( XAR-XBR ) 

GO TO 41 

40 YR2*Y-MX-XRI) /(TAN(THXAM) ) 

41 CONTINUE 
IF(YR2-5.5)42*42*43 

43 YR2»5*5 

42 CONTINUE 

DO 29 JJ»1*44 
CJ = JJ 

YRJ=(CJ-.5)*#25-5*5 

1F(YRJ-YL1)45»44*44 

44 IF(YR2-YRJ)45*46*46 
46 H«HIN 

XDS*HX-XRI )/H)**2 
YDS=((Y-YRJ)/H)**2 
PREF=GA0»GAM*RHM*RH0*DX2*DY2/PI 
DEN»(RHM + XDS-t-YDS )*#2 
REFL«PREF*DDlFni*JJ)/OEN 
IF(IWARN-l) 450*460*460 
460 REFL*2.*REFL 
450 DIR{ I *J)»OIR{I*J)-REFL 

45 CONTINUE 
29 CONTINUE 

36 CONTINUE 
28 CONTINUE 

WRITE(6*610)D1R( I *J) 

610 FORMAT (20X*‘ REDUCED INCIDENT FLUX IS '.E16.8I 
104 CONTINUE 
103 CONTINUE 
STOP 
END 

SUBROUTINE F I CT < X *Y *XL *6A*0n 

PI»3*1415926 

PSQ=1*/(PI**2) 

G=( (X-XU**2 + 1*H»»«5 
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ZEP*GA»( Y+.5) 

ZEM=GA»< Y-#5J 
6S»G**2 
ZPS«ZEP*»2 
ZM5=ZEM*»2 
PRE=PSQ/ i 2.*G) 

T1»ZEP/(GS+2PS) 

T2— ZEM/ (GS+ZMS) 

T3»ATANIZEP/G)/G 
T4— ATAN(ZEM/G)/G 
DI »PRE<M T l+T2-«-T3-*-T4 ) 

RETURN 

END 

SUBROUTINE FAST ( X lY »XLAfY 1 *Y2 tH*DLP ) 

PI*3*IA15926 

PSQ=l*/(PI**2) 

XSI»(X-XLA»/H 

G»(XSI**2+1*)***5 

Z1=Y-Y1 

Z2»Y-Y2 

GS=(H*G)**2 

ZS1»Z1*»2 

ZS2»Z2**2 

PRE=-PSQ/( 2.*G) 

T1»H*X2/IGS+ZS2) 

T2«-H*Z1/ {GS+ZSl ) 

ARG3=22/(H*G) 

ARG4»Z1/(H*G) 

T3=ATAN(ARG3)/G 
T4— ATAN( ARG4»/G 
DLP»PRE<MT1-H2+T3+T4) 

RETURN 

END 

//GO.SYSIN DO * 
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8.875 

9.375 
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DESIGN OF A QUARTZ LAMP THAT PRODUCES 
UNIFORM DIRECT RADIATION 


It would be desirable to develop a lamp geometry that provides a uniform 
incident flux on the surface beneath it. A system of that type would not be 
entirely possible if attention is restricted to cylindrical lamps which radiate 
symmetrically. However, it may be useful to determine the intensity distribu- 
tion leaving a lamp of finite length, which will provide the most nearly uni- 
form intensity distribution directly beneath it. The intensity distribution 
leaving the lamp could be controlled by varying the number of turns per centi- 
meter in the tungsten filament. That is, if the tungsten wire is of constant 
thickness, the intensity leaving the lamp will be direcly proportional to the 
resistance per centimeter, which is directly proportional to the number of 
turns per centimeter. 


From equation (2), page 4, the intensity distribution directly beneath 
the lamp is given by 

h Ii(i) 

. ~ ' D " r / 

0 


T, X i2(o,y) H . _ 
' D - 

n 


-4 [1 + - lyy 


<1^ 


(Bl) 


If Ijj is selected so that I(y) = r is the desired intensity distribution 
beneath the lamp, the problem would be reduced to finding Ii(c) so that 

(B2) 


1 = 


ih 


‘Inr 


5^ C 


[1 + Y^(c - IVV 

£ 

Equations of this type are known as Fredholm integral equations of the first 
kind (ref. 2). This particular problem could be solved by assuming that 
takes the form 

Ii(4) = ao + + 32?^ + 33^3 + . . . (B3) 

Under this assumption, equation (Bl) could be rewritten 


1 


u -n 


n^"-i [1 + 


di 


(P4) 


or 




(B5) 
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where 


' Xd * /u - 

Integrals in the form of equation (B6) could be integrated and yield 




(B6) 


h (^) = ^ ISjLil) — 


_ U - h) 

)2 1 + - %) 


^ tan~^ y(^ + h) 


tan 


y(c - h) 

Y 


(B7) 


(B8) 


(B9) 


u / N - 1 < 1 + + y _ Li-xlLerJik 

^ ' 2y2|i + y2(^ + h)^ 1 + y"(c - 

+ Ycftan'^ y(? + h) - tan"^ rU - ^)]j 

u / ^ - 1 ^ ^ ~ 

^ ' 2y^( 1 + ^ 

+ ]Li-l3S-i|f;an~^ y(? + %) - tan ^ y(^ ~ *5)]| 

and so on. 

The obvious problem which remains to be solved is the calculation of 
coefficients in equation (B4). Morse and Feshback (ref. 3) have shown 

that the series"in equation (B3) is not a useful form because of the difficulty 
in evaluating the a^. Rather, the series should have been written in the form 


of 

ii(c) = ao + ai(; + aio) + a2(^2 + 0121 C + ^ 20 ) 
+ a3(^^ + 01324^ «3 i 4 + “3o) + • • • 

Under these conditions, equation (B6) takes the form 
1 = aoho(c) + ai + “10^0(4)] ^ 


(BIO) 


(BID 


V W - L “^1 

+ a2rh2(?) ■*■ “21^1(4) ■*■ a2oi^o(^)J + • • ■ 

The basis for this polynomial expansion is that the Gram-Schmidt orthogonal iza- 
tion procedure could be used to construct orthogonal functions i|i^(c) f^om the 


h^(c) by 


<l'o(4) = ho(^) 

= hi(^) + aioho(4) 

'{'2(4) ~ ^ 2 (^) + + 0120*^0(4) 
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n-1 



If the inner product is defined by 

(f.g) = /^f(i)g(c)dc, 

^o(^) ^i(c) 3^® orthogonal if 

° " (hi.i^o) + «ioikoll^ 

-k 


(B12) 


(B13) 


(B14) 


where 

ll^oP “ /^'Po^(4)dE 
-h 

Obviously, equation (B14) would be satisfied if aio were given by 
~(hi >4 ^q) 

' UoF 

Similarly, 021 and a 2 o s*'® given by 

-(h2.’h) 


(B15) 


(B16) 


“ ll^oP 

These quantities were evaluated by numerically integrating 
inner products. 


(B17) 


the appropriate 


The magnitudes of ni„ and an are negligibly small when compared with 
the numerical integration accuracy. The value of 020 It a function of y 
and is tabulated below: 
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OPTIMIZATION COEFFICIENTS FOR UNIFORM DIRECT RADIATION 


Y 

“20 

^0 

2 

-.070536 

1.769 

4 

-.067466 

2.947 

6 

-.068425 

4.186 

3 

-.069984 

5.433 

10 

-.071436 

6.706 

20 

-.075908 

13.054 

40 

-.079131 

25.777 

60 

-.080364 

38.507 

80 

-.080994 

51.239 

100 

-.081360 

63.975 


32 

"0 

B 

13.100 

0.491 

36.864 

11.704 

2.157 

5.426 

10.257 

3.484 

2.944 

9.702 

4.764 

2.037 

9.434 

6.032 

1.564 

9.092 

12.364 

0.7354 

9.058 

25.060 

0.3615 

9.086 

37.777 

0.2405 

9.148 

50.498 

0.1312 

9.290 

63.219 

0.1469 


If the uniform intensity distribution is to be approximated by 

ii(c) ^ ao 

it is necessary to regroup the first three terms in equation (BIO) 
if ii(c) is approximated by 

ii(c) ^ ao + 92(0120 + 

where ao, aj, and a 2 are given by 


ao = 


ai 


and 32 = 


('l'2 ) 


(B13) 
That is, 

(B19) 

(B20) 


an approximate set of values for ai, and ^2 could be developed. Numer- 
ical integration of the expressions in equation (B20) shows that ai is neg- 
ligibly small and that ao and a 2 vary with y* These quantities along with 


ao 


are also tabulated above. (Note that ^2 = a 2 -) 

If ii(?) is written 
Ii(?) = ao(l + 3^^) 

the relative importance of the quadratic term is indicated by B and is shown 
in the table above. It can be seen that for large values of y (long lamps 
close to the surface), the nonuniform term becomes small; whereas for small 
values of y (short lamps away from the surface), the uniform contribution 
is small and the lamps should essentially be unheated at their centers. 
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